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Abstract. Despite its mathematical complexity, the mul- 
tiple gravitational lens can be studied in detail in every 
situation where a perturbative approach is possible. In 
this paper, we examine the caustics of a system with a 
lens very far from the others with respect to their Ein- 
stein radii, and a system where mutual distances between 
lenses are small compared to the Einstein radius of the to- 
tal mass. Finally we review the case of a planetary system 
adding some new information (area of caustics, duality 
and higher order terms). 

Key words: Gravitational lensing; Binary stars; Plane- 
tary systems 
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1. Introduction 

The investigation of a gravitational lens composed by a 
discrete set of point-like masses was opened by Schneider 
& Weifi in 1986. They considered a binary lens, deriv- 
ing analytical expressions for caustics in the case of lenses 
with the same mass. Erdl & Schneider (1993) extended 
their analysis to an arbitrary mass ratio, confirming the 
presence of three possible topologies in the caustic struc- 
ture. 

Besides these two fundamental papers, many other 
contributions have been given to this field by several au- 
thors studying particular aspects of this matter. It has of- 
ten been considered the relation between the binary lens 
and the Chang-Refsdal (1979; 1984) model describing the 
behaviour of a point-like lens immersed in a uniform grav- 
itational field (Schneider & Weil3 1986; Dominik 1999). In 
the special case of very asymmetric lens, as that of a planet 
orbiting a star, the Chang-Refsdal approximation can give 
very reliable results (Gould & Loeb 1992; Gaudi & Gould 
1997). In the search for good approximate models, also the 
quadrupole lens was given a role in close binary systems 
(Dominik 1999). 

However, while the picture for the binary lens is some- 
what well defined at present day, the same cannot be said 



about general multiple lenses. The attempts to enlighten 
this kind of lenses have really been very few. This is be- 
cause of the mathematical complexity of the lens equation. 
Yet, multiple lenses can play an important role in several 
cases and therefore are not just mathematical curiosities. 
Planetary systems are likely to be composed by a star 
with more than one planet, so it is important to study 
them in their full complexity in order to interpret plan- 
etary microlensing events correctly, especially those with 
high magnification (Wambsganss 1997; Gaudi, Naber & 
Sackett 1998; Bozza 1999). N-body systems made of sim- 
ilar masses are hard to be found as gravitationally stable 
systems, yet the projected distances between far stable 
subsystems could be compatible with the appearance of 
"interaction" in the gravitational lensing behaviour; some 
considerations about microlensing light curves by multiple 
lenses have been done by Rhie (1996). 

Another context where multiple lensing is surely of 
great interest is cosmological lensing by rich galactic 
clusters: as a first approximation, all lensing galaxies 
can be considered as point masses, then the multiple 
Schwarzschild lens constitutes a good starting point for 
the analysis of caustics and images produced by such sys- 
tems. 

The aim of this work is to give a systematic inves- 
tigation of the caustics of multiple lenses whenever they 
can be referred to well-defined situations to be employed 
as starting points in perturbative expansions. The three 
cases where this happens have been pointed out by Do- 
minik (1999) and are resumed here. 

In Sect. 2 the lens equation for multiple lenses and the 
basic definitions are given; Sect. 3 deals with the case of 
a lensing mass that is very far from the others; in Sect. 
4 a system formed by near (with respect to the Einstein 
radius of the total mass) lenses is studied; Sect. 5 recon- 
siders the case of a multiple planetary system, already 
examined in a former work (Bozza 1999); Sect. 6 contains 
the conclusions. 



E-mail valboz@physics.unisa.it 



Valeric Bozza: Caustics in special multiple lenses 



2. Basics of multiple lensing 

Let 
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be the Einstein radius of a mass Mq placed at distance 
DoL from the observer and at distance Z3ls from a source 
having distance Z?os = -Dql + £'ls from the observer. 
The coordinates in the lens plane normalized to R^ will 
be denoted by x = (xi;a;2), while the coordinates in the 
source plane normalized to -Re'H^^ ^^"^ y = (j/i;2/2)- AH 
masses are meant to be measured in units of Mq which 
is a typical mass of the problem we are considering. For 
example it can be the solar mass if the lensing objects 
are stars, or a typical galaxy mass if we are considering 
extragalactic lensing. 

According to the standard theory of gravitational lens- 
ing, the deviation induced by a system formed by n masses 
mi, ..., m„ whose projections in the lens plane are at po- 
sitions xi, ..., x„ is the sum of the deviations produced by 
the single objects. So the lens equation reads: 
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Given a source position y, the x's solving this equation 
are called images. Their properties can be studied through 
the Jacobian matrix of the lens application. In particular, 
the determinant of this matrix is: 
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where A^ = (A^; Ai2) = x-x^. 

The magnification of an image is given by the absolute 
value of the inverse of det J calculated at the image po- 
sition. Whenever the Jacobian determinant vanishes, an 
image at that point would be (in ray optics) infinitely am- 
plified. These points are thus said to be critical and their 
corresponding points in the source plane, found by use of 
Eq. (Q), constitute the caustics. 

Caustics appear as closed cusped curves and split the 
source plane in several domains. The number of images 
corresponding to a single source only depends on the do- 
main where the source is located. In neighboring domains 
the number of images always differs by two units. The 
caustics of a given lens model provide very useful hints for 
the description of the general properties of the system so 
that their study is essential to understand the behaviour 
of a lens. 

To find the critical curves and then the caustics of a 
multiple lens, one has to solve the equation det J — 0, 
which is very complicated. However, if it is possible to 



single out some parameters which are very large or small 
with respect to the others, a perturbative approach can 
be tried (Bozza 1999). The situations we are to explore in 
this paper always refer to the single Schwarzschild lens as 
the zero order solution. The critical curve of a lens with 
mass m placed at the origin is a circle with radius \/m, 
while the caustic is reduced to the origin itself: y = 0. 

3. Single lens perturbed by far masses 

Let's start considering an isolated point-lens. The presence 
of other masses that are very far from the first one (with 
respect to all Einstein radii) warps the circular critical 
curve, giving rise to an extended caustic. 

Let's consider a mass mi placed at the origin of 
the lens plane. The other masses 7712 , ■ • ■ , Wn are at the 
positions X2,...,x„. In polar coordinates, we put x^ — 
(pj cos If j^; Pj^ sinyjj) and the generic coordinate in the lens 
plane is x = (r cos 1?; r sin d) . The hypothesis we start from 
is that pj 3> y/rnj for each i and j. This allows us to 
expand the Jacobian determinant (0) in series of inverse 
powers of p^. If we stop at the fourth order we have: 
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We see that a linear superposition principle is valid at 
the second and the third order. In the fourth order there 
are "interaction" terms that make the calculations more 
difficult. 

The zero order solution is simply the Einstein radius 
r = ^Jm\. We build the general solution as an infinite 
series in Xj p{. 



£2 +£3 + ■••) 
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where Sj ~ l/p^- Substituting this expression in Eq. (|4|) 
and expanding again, we can easily solve the equation of 
critical curves at each order to find the perturbations. The 
first three orders are: 
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Recalling that the shear produced by one mass in a 
Chang-Refsdal approximation is 7^ = ^, the second or- 
der can be easily recognized as the sum of the Chang- 
Refsdal expansions for each mass to the first order in the 
7j. Actually, the full second order can be identified with 
a Chang-Refsdal lens expanded to the first order in this 
total shear: 



n 



=2 



rrii sin {2ip^ 



y—^ rui cos(2(^j) 

,1=2 



Pi 



(7) 



and oriented along the direction at the angle: 
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So, the Chang-Refsdal lens provides a first approxi- 
mation for the critical curves of a far lens not only in 
the binary case (Dominik 1999) but even for an arbitrary 
number of lenses, as already stated in (Chang & Refsdal 
1979). 

Now we exploit the perturbative results just shown to 
find the caustic formed in this case. It suffices to apply the 
lens equation (|2|) and expand it to the third order again; 
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The first term is independent of i? and only fixes the 
position of the caustic which is displaced towards the other 
masses; the two second order terms describe the shape of 
the caustic reproducing the Chang-Refsdal limit; the third 
order terms correct this approximation. 

The Chang-Refsdal model describes the gravitational 
lensing by a mass embedded in a uniform gravitational 
field. For 7 > 1, the caustics of this model are diamond- 
shaped curves with four cusps which are independent on 




Fig. 1. The caustic of a lens perturbed by two other bod- 
ies. The masses are mi = 0.25, m2 = 0.25, rn^ — 0.5. The 
positions of 7712 and 7713 are P2 = 4:, (p2 = and p^ = v20, 
(/Sg = arctan 0.5 respectively. The solid line is the pertur- 
bative caustic to be compared with the dashed line that is 
the exact one, almost completely hidden behind the solid 
line. 



the versus of the gravitational field. The approximation of 
uniform field is reasonable in our case. However, it breaks 
down quite soon when the distances are not so high. It 
happens that our lens begins to feel the non uniformity of 
the field and the caustics lose their symmetry, elongating 
in the direction of the other masses. The third order de- 
scribes this elongation very well, considerably extending 
the range of applicability of the perturbative results. 

The breakdown of this perturbative expansion rises 
when the Einstein ring of the first mass is too close to 
the critical curves produced by the other masses. A par- 
ticular case is obtained when several small masses are very 
close each other. The critical curve generated by such a 
subsystem has a radius of the order of the square root of 
its total mass. So even if the condition p^ ^ J"^] is sat- 
isfied for each single mass, the distance of the first mass 
to the subsystem could be smaller than the radius of its 
total critical curve, causing the failure of the perturbative 
hypothesis. 

In Fig. ^, we see an example of the power of the per- 
turbative approach in reproducing the numerically found 
exact results with great accuracy not taking care of the 
number of masses. 

Fig. shows an interesting comparison among the dif- 
ferent orders approximations and the exact caustic for 
different distances between the two lenses. We see that 
even at separations about 10-\/2 times the Einstein ra- 
dius of each lens, the second order approximation (Chang- 
Refsdal) seems quite poor, while the third order caustic 
is practically coincident with the exact one up to five 
Einstein radii. As previously said the inadequacy of the 
Chang-Refsdal approximation is in its impossibility in giv- 
ing account of the asymmetry of the caustic, which is al- 
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Fig. 2. Caustic of a lens with a far companion of the same mass (0.5), along the direction ip^ = 0- The dotted curves 
are the exact caustics, the dashed curves are the second order perturbative ones (Chang-Refsdal) and the continuous 
lines are the third order ones, (a) p^ — 10, (b) p,^ = '''■5, (c) P2 = 5, (d) p^ = 2.5. 



ready evident at quite far separations. This justifies us in 
considering the third order that is the lowest order repro- 
ducing this asymmetry. 

Now let's turn to the cusps of these caustics. These 
points are characterized by the vanishing of the tangent 
vector: 



^2 (^) = 



(10) 



At the second order we have a Chang-Refsdal caustic 
rotated at an angle <~p given by Eq. (ph from the yi-axis. 
Without loss of generality we can put this angle to zero. 
Then it is well known (Dominik 1999) that the curve is 
symmetric for reflections on two axes passing through the 
center of the caustic and has four cusps at the intersections 
with these axes, corresponding to z9 = 0, 5, tt, ^. 

At the third order, summing and subtracting the two 
Eqs. (p^), we get: 
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The two Eqs. can be simultaneously satisfied only if 
F {•&) vanishes. For a binary system, this equation has six 
roots: 

i? = 
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The first two are the cusps along the j/i-axis already 
present at the second order. The second two are the mod- 
ification of the other two Chang-Refsdal cusps (in fact, 
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The result is: 




Fig. 3. A particular of Fig. ||d showing the two additional 
cusps yielded by the perturbative expansion. 

when P2 tends to infinity, they approach ±-|). The last 
two cusps are imaginary for pj > 4^r7ii, but become real 
for lower distances, giving rise to a "butterfly" geome- 
try around -d = tt (Fig. ||). For binary systems, this is 
not what happens in exact results which always yield four 
cusps; yet, with more than two lenses, an increment in the 
number of cusps is effectively present when the separation 
between the lenses is not very high. 

Having at our disposal good approximate expressions 
for the caustics of multiple lenses, an interesting quan- 
tity to compute is the area covered by these curves in the 
source plane. We do this calculation for an arbitrary mul- 
tiple lens at the second order and for a binary system at 
the third order. 

For the Chang-Refsdal approximation, again we as- 
sume that the angle ip defined in Eq. (^) is zero. Then the 
caustic is oriented along the yi-axis and its area is given 
by: 



Vl 
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We have previously expressed j/2 and yi as functions 
of ■&. Considering that, when ^ = 0, j/i is y'^^^ and, when 
t9 = TT, j/i is y""", also considering that 2/2 is negative for 
< i9 < TT and positive for tt < i? < 27r, the area can be 
calculated by the following integral: 
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Remembering the expression (0) for 7 in the physical 
system we are analyzing, we see that the main depen- 
dences of the area of the caustic on its parameters are 
quadratic in the masses and inverse fourth power type in 
the distances. 

For the third order calculation, a direction of symme- 
try can no longer be determined unless we limit ourselves 
to the binary system. The procedure just exposed can be 
followed for p^ > Ay^mi, when the geometry of the caus- 
tic remains unaltered. Then, we obtain a correction to the 
previous result: 



AA = A-K- 



P\ 



(17) 



So, not only the third order stretches the caustic but also 
adds a positive contribution to its area. 



4. Close multiple lenses 

Now we consider a set of n point-like lenses whose mu- 
tual distances are small compared to the total Einstein 
radius of the system. The starting point for our analy- 
sis is the Schwarzschild lens we would obtain if the total 
mass M were concentrated at the barycentre. The sepa- 
rations between the masses modulate the deviations from 
the Schwarzschild lens and then constitute the perturba- 
tive parameters in our expansion. 

Now there is no privileged mass, so the most reasonable 
choice of the coordinates origin would be the centre of 
mass. We derive our results in the general case but we 
will choose the centre of mass system for some particular 
considerations. 

Besides the main critical curve, which is a slight mod- 
ification of the Einstein ring of the barycentral lens, some 
secondary critical curves are present near the centre of 
mass. They give rise to very small and far caustics which 
must be treated separately. 

^.1. Main critical curve and central caustic 

We consider our masses ?7ii,...,m„ placed at the po- 
sitions xi,...,x„. In polar coordinates, we put x^ = 
(/9j cos(y9j; Pj siniy9j) and the generic coordinate in the lens 
plane is x = (r cosi?; rsin??). Now we assume that p, <^ 
\fM for each i. We carry the expansion of the Eq. (||) in 
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series of powers of p^ up to the second order: 
det J = 1 J — > PiTTii cos [it — ip^j + 
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3 ^ m.Mpl cos (2?9 - 2(^J + 



2M 

^•6 



+2 ^ m^nijp^pj cos ((^^ - (^j) 



(18) 



The zero order terms give the Einstein ring for a point 
mass M at the origin. We see that interference terms arise 
at the second order and this is a good reason to stop our 
expansion to avoid long calculations at higher orders. 

We can consider an expression analogue to Eq. (||) for 
the radial coordinate r: 



/M(l + ei +S2 + -) 
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where now Ej ~ p^ . Substituting in Eq. dlq) and expand- 
ing again, we can solve for the Ej at each order: 

n 

£1 = j^ E "^iP* COS {■& - ip.^) 

£2 = 4T73 1 6 f: Mm^pf COS (2i9 -2ip,) + 



" E rrumjPiPj [5 cos (2i9 - cp, - ^Pj) + cos {^p^ - ip^)] \ 

ij = l J 

(20) 

Now that the perturbations to the critical curve are 
known (to the second order), we can find the caustic by 
use of the lens equation: 
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The first order perturbations only displace the caustic 
to the centre of mass. Choosing this point as the origin 
of the coordinate system, we can eliminate these terms. 
The second order perturbations then yield the first non 
trivial deviations from the Schwarzschild lens. These per- 
turbations are in the form of products between the masses 
forming the set of lenses, so there is no superposition prin- 
ciple here. 



This perturbative expansion is tightly related to the 
quadrupole lens, as pointed out by Dominik (1999) for the 
binary lens. The quadrupole lens equation has the form: 



where 
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is the quadrupole moment. 

A careful analysis of the second order Jacobian of the 
close multiple system (18) in the centre of mass system 
reveals its coincidence with the Jacobian of the quadrupole 
lens expanded to the first order in the quadrupole moment, 
in polar coordinates. The elements of the matrix (G3) are: 



Qi = M J2 iTT'tPi cos {2ipi) 

4=1 

n 

Q2 = M J2 "m^pf sin (2<y9,) 

i=l 
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The two eigenvalues of the quadrupole matrix have the 
same absolute value 



Q^JQ\ + Ql 



(25) 



and opposite signs. We can define the orientation of the 
lens as the direction of the eigenvector corresponding to 
the positive eigenvalue, because it reduces to the positive 
xi direction in the case of a binary system with the masses 
aligned on the xi-axis. The angle that this eigenvector 
forms with the a: 1 -axis is: 



If = arctan 



Q-Qi 
Q2 



(26) 



The caustic of a quadrupole lens is a diamond-shaped 
figure symmetric for all reflections on the eigenvalues axes. 
If we choose the coordinates so as to have the caustic ori- 
ented along the xi-axis, the four cusps are at the positions 
■d — 0, ^,7r, ^. Fig. ^ shows some examples of central 
caustics for a binary system with a mass ratio q = 3. The 
second order approximation works well up to separations 
of tenths of the total Einstein radius. It is interesting to 
remark that the second order result is always symmetric 
even if the system of lenses has no symmetry at all. The 
symmetry is only lost when the perturbative hypothesis 
begins to be more forced, as we see in Fig. Ha. 

For the extension of the central caustic of a system of 
close multiple lenses, we can directly refer to a well ori- 
ented quadrupole lens caustic expanded to the first order 
in Q. The area of such a caustic can be calculated with 
the same procedure explained in the previous section for 
Chang-Refsdal caustics. The result is: 



A = 



37rg2 
2M3 



(27) 
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Recalling Eqs. (|2J) and (p5|), we see that the area goes 
as the fourth power of the distances between the masses, 
while the mass dependence is not so immediate because of 
the presence of the cube of the total mass in the denomi- 
nator. 

4-2. Secondary caustics 

In the neighbourhood of the centre of mass of the system, 
some critical curves in the form of small ovals can appear 
(see Fig. ||). The number of these curves is not fixed by 
the number of lenses but depends on the position and the 
values of the masses composing the system. The binary 
lens is the only case where the number of critical curves 
is always fixed to two. 

The perturbative approach can be used in the case of 
the binary lens to find compact formulae for the positions 
of the secondary critical curves and consequently the po- 
sitions of the caustics. The general multiple lens involves 
high degree equations which do not allow any analytical 
considerations. 

Let's consider two masses mi and 7712 separated by a 
distance a along the xi-axis. In the centre of mass system, 
the equation det J = can be rationalized in the following 
form: 
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where M = mi -\- mi- 

We search for solutions of the first order in the sepa- 
ration a. So we put xi — aAxi, X2 — aAx2 and expand in 
powers of a. The first non trivial order is the fourth, so the 
equation is of the fourth order in Axi and Aa;2. Introduc- 
ing the polar coordinates Axi — Ar cos f?, Aa;2 = Ar sin "(9, 
our equation can be solved for Ar. The four solutions 
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Ar = ± — [{m,2 — mi) cos-i? + \Jmi1m2 sin^- 



±i (Y^mim2C0S?? + (mi — m2)sini?)] (29) 
are real only if 
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This equation is satisfied by the angles: 
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1? = ± arccos ■ 
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Inserting these values in Eq. (M), discarding the negative 
solutions, and returning to the cartesian coordinates, we 
have: 
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If we put these values in the original equation, there are 
only two acceptable solutions: 
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These two points represent the centre of the small ovals at 
the first order in the separation between the two masses. 
The shape of the curves can be studied only at higher 
orders but high degree polynomials heavily complicate the 
calculations. 

Through the lens equation, we can find the positions 
of the corresponding caustics: 



yi = 



2/2 = ± 



TOi — m2 



a 



2^mim,2 



(34a) 
(34b) 



We notice that when a tends to zero, the two caustics 
become infinitely far. 

These formulae are compatible with those for the posi- 
tions of the couple of planetary caustics of planets internal 
to the Einstein ring, given in (Bozza 1999), in the limit of 
planets very close to the star. 

5. Planetary systems 

A particularly interesting case of multiple lensing is that 
of planetary systems. Here a stellar mass is surrounded by 
planets having masses much smaller. Their effects can usu- 
ally be treated as slight perturbations to the main lensing 
object. 

The structure of the caustics of this lens has been 
widely explored by perturbative methods in a previous 
work (Bozza 1999). Here we shall recall the main formulae 
to complete the picture of perturbative results in multiple 
lensing and make some additional considerations. 

We can distinguish between the central caustic which 
is generated by the distortion of the star's Einstein ring 
and the planetary caustics which are the images of the 
planetary critical curves and can be well approximated by 
Chang-Refsdal critical curves. 

5.1. Central caustic 



rriimi 



(31) The central caustic of a star with mass mi placed at 
the origin, surrounded by planets with masses TO2, ■•■,m„ 
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placed at X2, ...,x„ respectively is expressed by the follow- 
ing parametric form: 



j/i (t?) = 2^/m:^e (d) cos I? - X^ 



miA°^ 



% [(AO0^ + (AO,)^] 



yi if) = 2^/^e {i}) sin I? - X] 



^, {Ao,r+iAo,r 



(35a) 



(35b) 



where 



e{i}) = -cos2i?^; 



(A?r) 



(^2)' 



[{AO^f + {A^,y 



3in 2z9 Y^ 



A" A° 



2 [(A«i)' + iAO,y 



(36) 



and A^ = (^mi cosz? — x^; y^mT sin i9 — a;i2)- 

In many studies of the central caustic of a planetary 
system, a principle of duality between planets external 
and internal to the Einstein ring (hereafter, we shall sim- 
ply refer to them as external and internal planets, respec- 
tively) was often claimed on the basis of the observation of 
the shape of numerical caustics (Griest & Safizadeh 1997; 
Dominik 1999). This principle can be directly verified on 
these formulae which are invariant under the transforma- 
tion 



TOi 



(37) 



So the conjecture of duality has an effective analytical 
basis. 

The central caustic is generally a self-intersecting 
curve, except for the case of the single planet, when it has 
the already encountered diamond shape. Then, its area 
can be calculated with the same method of the previous 
sections. The result is: 



A 



4toi 



1 



(mi + pI) {m,l - Amipl + p|) 



mi 



2|3 

P2I 



(38) 



It positively diverges when P2 -^ \Jm\. This hap- 
pens because in this limit the perturbative approach is no 
longer valid and the fusion between the central and the 
planetary caustic occurs. The area vanishes when P2 ^ 
or y02 ^ 00 because, in these limits, the Schwarzschild lens 
is recovered and the central caustic reduces to a point. 

5.2. Planetary caustics 

The planetary critical curve is always localized in the 
neighbourhood of the planet and assumes the shape of 
an elongated ring, when the planet is outside of the star's 
Einstein ring, or splits into two specular ovals, when the 
planet is inside. 

The zero point of the expansion of these critical curves 
is then the position X2 of the planet we are considering. 



The first non trivial order is x — X2 ^ \A«^, since the 
critical curve of a very far planet is nothing but its Einstein 
ring with radius ^m^- 

In (Bozza 1999), the first order planetary critical 
curves were derived and their relations with the Chang- 
Refsdal ones was discussed. It is interesting to go farther 
in the expansion to understand how this relation is broken 
and to see the effects of other planets. 

It is useful to consider polar coordinates around the 
planet position X2. Let's then set x— X2 — (r cos ^; r sin-iJ). 
Now we expand the radial coordinate r: 



£1 +£2 



(39) 



where e,- 



i/2 

-i ' 



The first two orders of the equation det J = give: 



2vn\mi cos (2-d — 2ip2 



P\ 



^\pI 



4m2 



e-i 



4:17111712 cos (2i9 — 2(^2 
SiPI 



-£2- 



4mi \mie\ cos {-d — 1^2) + "^-2^2 cos (3?? — 3(^2)] 



£lP2 



(40) 



We see that the first row, representing the lowest order 
Jacobian, is just the Jacobian of the Chang-Refsdal lens 
1712 with shear 7 = ^. Then the lowest order critical 
curve is exactly Chang-Refsdal: 



£1 



\ 



m2 



^ cos (2^ ■ 



■ 2^2)±\n 



■ sm 



(2^ - 2^2 



1- 



P2 



(41) 



According to the double sign, two branches are present. 
For external planets only the higher is real, while for inter- 
nal planets both branches are real in two intervals centred 
on 7? = ±f (Bozza 1999). 

We can notice that if we perform an expansion of Eq. 
(ETl) to the first order in the shear, we obviously get the 
same result of Sect. 3, with the roles of m2 and mi in- 
terchanged. However, in this section we are analysing the 
perturbations in the masses of the planets, so this expan- 
sion is not interesting for our purposes, though it gives the 
connection with the previous calculation. 

The full Eq. (^) can now be employed to find the 
second order perturbation: 



£2 = -- 



miE^ [mi£i cos {■& — (^2) + ™2P2 cos (3i? — •^^2)] 



m2p\ [mie\ cos (2t? — 2ip2) + "^2^2! 



(42) 
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As usual, through the lens equation we can write down 
the formulae for the caustics: 



2/1 W= (P2-^ cos^2+fei-ff cosi^ 



P2 



e2Cosz9 + 



~mi -3 

f2 



E™ P2 COSip2~Pi COSiPi 



i=3 



' P2+Pz'^P2Pi COs{ip2-iPi) 



(43a) 



y2 W = 



P2- '^jsinv32+ [^1 - ^)sin^ + 

misi sin(T9 — 2(^2) 



-nil 



P2 

E2P2 sin(' 



fl + ^jeasini?- 

-2<P2)-Ei sin(2i9-3y2) 1 



W 



E™ P2siny2-PiSiny. 



i=3 



' P2+Pi -2p2Pi cos(¥>2-Vi) 



(43b) 



At the second order a sum containing the effects of the 
other planets appears. It does not depend on d and then 
only represents a displacement of the caustic towards the 
other planets. 

Fig. P shows a comparison between the perturbative 
results and the exact curve for an external planet. The 
loss of symmetry of the figure, not present in the first 
approximation, is reliably followed by the second order 
curve. Moreover, it is to notice that the displacement of 
the caustic, due to other planets, is correctly reproduced. 

The critical curves of internal planets are not defined 
for all values of -d. This fact brings in divergences at the 
border points where the solution switch from imaginary to 
real. So the second order results are very good everywhere 
except for a neighbourhood of these border values of the 
angle. 

As for the other caustics, we can easily calculate the 
area of planetary caustics exploiting their symmetries. The 
case of external planets is analogue to the previous ones, 
so we can directly give the result: 



A = 2m2 
where: 



^ / TT 771? \ „ /tT 777? 



Fiip-m) = /[I- 



777 sm 



^] 



-1/2 



d-ff 



E{(p;m) 



I [l~m sin^ i9] ' 



dd 



(44) 



(45) 



(46) 



are the elliptic integrals of the first and the second kind 
respectively. 

The calculation of the area for the caustics of the inter- 
nal planets is not very different. We have to integrate on 
the higher branch and subtract the integral of the lower so 



as to obtain the area included between the two branches. 
The formula (Hq) must then be modified this way: 



(47) 



where the two branches are real. 
The result is: 



A 



1712 [2E[ 279^i„; ^ 



- 2E ( 27?„,ax; ^ 

P2 



-F ( 27?„,i„; 



F\ 2i9 



maxi T" 
P2 



(48) 



Of course, as there are two planetary caustics for internal 
planets, this number must be doubled. 

The area of the planetary caustics is of order 7772. In 
fact it is well known that they are much greater than the 
central caustic, which is instead of order rn^ (see Eq. (p8[)) 
(Dominik 1999). 

In Fig. we see a plot of the area of the planetary 
caustic as a function of the separation of the planet from 
the star. We see the divergence in correspondence of the 
"resonant lensing" situation (pj = \A^) where the per- 
turbative theory cannot be applied. 

6. Conclusions 

The results presented in this work prove the great fiexibil- 
ity of perturbative methods in the resolution of problems 
in gravitational lensing. By different choices of perturba- 
tive expansions, it is possible to obtain very useful analyti- 
cal approximations, which can be employed for further ap- 
plications. In the case of a discrete set of point-like lenses, 
we gain a deep understanding of the distortion effects with 
respect to the Schwarzschild lens. This is the only ana- 
lytical knowledge we have about the shape of caustics in 
multiple lenses. 

Carrying the expansions to higher orders allows 
some interesting considerations about the validity of 
some approximations such as the Chang-Refsdal or the 
quadrupole lens and the limits of applicability of a su- 
perposition principle for the effects coming from different 
masses. 

The calculation of the area of the caustics is another 
important possibility offered by the achievement of these 
analytical approximations. 
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Fig. 4. Central caustic of a binary lens with m\ = 0.75 
and 1712 = 0.25 and separation (a) 0.2, (b) 0.1, (c) 0.05. 
The dashed curve is the exact caustic and the solid is the 
perturbative one. 
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Fig. 7. Area of the planetary caustics as a function of the 
distance of the planet from the star. 
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Fig. 5. Secondary critical curves generated by four equal 
masses, represented by the small crosses in the figure. 
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Fig. 6. Planetary caustic of a jovian planet (m,2 = 10^"^) 
at position (1.2; 0). The caustic is perturbed by the pres- 
ence of another jovian planet at (0; 1.5). The dashed curve 
is the first order approximation, the solid one is the second 
order approximation, which is almost coincident with the 
exact one (dotted curve). 



